FCC-BCC-Fluid triple point for model pair interactions with variable softness 
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It is demonstrated that the coordinates of the fcc-bcc-fiuid triple point of various model systems 
are located in a relatively narrow region, when expressed in terms of the two proper variables, 
characterizing the softness and strength of the interaction force at the mean interparticle separation. 
This can be regarded as a consequence of the "corresponding states principle" for strongly interacting 
particle systems we have put forward recently [S. A. Khrapak, M. Chaudhuri, and G. E. Morfill, J. 
Chem. Phys. 134, 241101 (2011)]. The related possibilities to predict the existence and approximate 
location of the fcc-bcc-fiuid triple point for a wide range of pair interactions with variable softness 
are illustrated. Relation of the obtained results to experimental studies of complex (dusty) plasmas 
are briefly discussed. 

PACS numbers: 64.60.Ej, 64.70.D-, 52.27.Lw 
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Many systems of strongly interacting particles in ther- 
modynamic equilibrium are known to form face-centered- 
cubic (fee) or body-centered-cubic (bec) lattices in the 
solid phase. Normally, for sufficiently soft interactions 
there is a tendency to form a bee solid. On the contrary, 
steep interactions favor the fee structure. When the in- 
teraction potential exhibits variable softness, a polymor- 
phic phase transition between fee and bec solids can be 
present and the fcc-bcc-fiuid triple point emerges on the 
respective phase diagram. 

Finding the location of fcc-bcc phase transition re- 
mains a difficult and resource demanding (computa- 
tional) task since usually a very tiny difference between 
the free energies of the corresponding phases has to be re- 
solved. In this Letter we discuss a very simple approach 
to approximately locate the fcc-bcc-fiuid triple point. The 
approach is based on analyzing the shape of the pair in- 
teraction potential only and does not require any kind of 
numerical simulations. 

The qualitative similarities in the phase behavior of 
model systems with soft repulsive interactions have been 
discussed previously from various perspectives. In partic- 
ular, Prestipino et al. [l|, Q discussed the common topol- 
ogy of the phase diagrams of the Gaussian core model 
(GCM), inverse-power- law (IPL), and Yukawa model po- 
tentials in the vicinity of the fcc-bcc-fiuid triple point. 
In an attempt to relate the phase behaviors of the GCM 
and Yukawa potentials to that of the IPL model, they re- 
quired the logarithmic derivatives of the three potentials 
to match, for separations close to the average nearest- 
neighbor distance. As a result, the effective exponents 
for the GCM and Yukawa interactions at the triple point 
were shown to lie not too far from the corresponding IPL 
value. This was a useful observation, although perhaps 
insufficient for direct practical applications. 

As we conjectured recently similarities in the be- 
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havior of two different systems of strongly interacting 
particles can be expected when the interparticle force 
and its first derivative (evaluated at the mean interparti- 
cle distance) in one system are equal to those in the other 
system. The simplistic arguments behind this conjecture 
are as follows: In a strongly coupled state the particles 
form a regular structure where large deviations of the 
interparticle separation from its average value are very 
seldom. Consequently, the state of the system should 
be virtually insensitive to the exact shape of the inter- 
action potential at short distances. If, in addition, the 
potential decays sufficiently fast for distances beyond the 
mean interparticle separation, its long-range asymptote 
is of little importance for the thermodynamics of the sys- 
tem [H . The behavior of the potential at the average in- 
terparticle separation plays a dominant role. Physically, 
the force and its first derivative are the two quantities 
which should mainly matter. Using the IPL potential as 
a reference system, we identified the two scaled variables 
which conveniently characterize the state of a strongly 
coupled system of particles interacting via the pair po- 
tential U(r). These are the generalized softness param- 
eter s = [-1 - C/"(A)A/J7'(A)] _1 and the reduced force 
T = —U'(A)A/T, where A is the mean interparticle sep- 
aration and T is the temperature (in energy units). The 
approximate "corresponding states principle" postulates 
that two different systems of strongly interacting parti- 
cles having the same values of s and T behave alike 0]. 
As a test of the vitality of this principle, we demonstrated 
that the melting curves of various model systems such as 
IPL, Yukawa, conventional Lennard- Jones (LJ), n — 6 
LJ, exp— 6 and (slightly less successfully) GCM are es- 
sentially collapsing on a single curve when plotted in the 
(s,J-) plane and, hence, the "universal melting curve" 
emerges Q. 

The question we address here is whether this corre- 
sponding states principle can also be useful to estimate 
(at least approximately) the location of the fcc-bcc-fiuid 
triple point. In order to answer it, we carefully ana- 
lyze the available data on the fcc-bcc phase transition 
for three selected repulsive potentials. We show that the 
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triple point coordinates in the (s, F) plane are indeed rel- 
atively close for all the systems considered. In addition, 
the phase diagrams in the vicinity of the triple point are 
topologically equivalent. 

We consider the IPL potential, U(r) = e(a/r) n , 
the Yukawa (exponentially screened Coulomb) poten- 
tial, U(r) = e(cr/r) cxp(— r/er), and the GCM potential, 
U(r) = ecxp(— r 2 /<7 2 ). Here e and a are the energy and 
length scales, and r is the distance between the two in- 
teracting particles. The IPL potential defines a system 
of "soft" spheres which can serve as a model for a wide 
range of physical systems, including the one-component 
plasma (OCP) and simple metals under extreme thermo- 
dynamic conditions. At n — > oo it approaches the hard 
sphere (HS) limit. The Yukawa (Dedbye-Hiickel) poten- 
tial is extensively used in the context of colloidal sus- 
pension, physics of plasmas, and complex (dusty) plas- 
mas [sMiol] . The GCM potential [HI is frequently used 
in the context of soft matter physics, e.g. to describe 
effective interactions in many body systems of polymer 
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The phase state of the considered systems in ther- 
modynamic equilibrium can be conveniently character- 
ized by the reduced density p* = pa 3 and temperature 
T* = T/e (For the IPL potential, the single variable 

p*T* fully describes the thermodynamic state of the 
system). Another possible choice of state variables for 
purely repulsive interactions uses the coupling parame- 
ter r = U(A)/T and the screening parameter k = A/er, 
where A = p -1 / 3 is the structure- independent interpar- 
ticle distance. This choice is particularly suitable for 
Yukawa interactions and is commonly used in the field 
of complex (dusty) plasmas [f| 0] ■ 

The phase portraits of the IPL, Yukawa, and GCM sys- 
tems are relatively well studied 0, [H, [Ur{27| ■ Concerning 
the fcc-bcc-fluid triple point, Agrawal and Koike [l6[ esti- 
mated its location for the IPL system as p*T* ~ 2.173 
(in the fluid phase) and 1/rt ~ 0.16. They noted, how- 
ever, that this estimate is subject to significant uncer- 
tainty. Later, Prestipino et al. [2[ performed more de- 
tailed study of the IPL phase diagram near the triple 
point. Their result is p„T* 3 ^ n ~ 1.787 (in the fluid 
phase) and 1/n ~ 0.14. In the same work the triple point 
temperature of the GCM system has been estimated as 
T» ~ 0.0031. The corresponding fluid phase density is 
~ 0.093. Dupont et al. [2l| reported the following 
triple point coordinate for the Yukawa model: k = 6.75 
and T* ~ 0.43 x 10~ 4 . Later, Hamaguchi et al. (HJ deter- 
mined the location of the triple-point in Yukawa systems 
as k = 6.90 and T ~ 3.50. Finally, in a recent paper by 
Hoy and Robbins [24| , the triple point of Yukawa systems 
has been located at k = 7.70 and T ~ 3.11. For consis- 
tency, we take the fluid density at the triple point where 
appropriate. For the data related to the Yukawa poten- 
tial, the triple-point density has a single value, defined 
by the estimated value of k. 

The corresponding values of Ttp , .Ftp , and stp at the 



TABLE I: The estimated values of the parameters T, T and s 
at the fcc-bcc-fluid triple point for various model interaction 
potentials. 



Potential 


Ttp 


Ftp 


s T p 


Ref. 


GCM 


2.47 


24.08 


0.129 


[2] 


IPL 


3.98 


28.46 


0.14 


m 


IPL 


5.05 


31.57 


0.16 


[16] 


Yukawa 


4.03 


31.26 


0.145 


[21] 


Yukawa 


3.50 


27.65 


0.142 


[22] 


Yukawa 


3.11 


27.09 


0.128 
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triple point are summarized in Table HI The scattering in 
the reported Ttp values is relatively large, up to a factor 
of two (see the first column in Table HJ. On the other 
hand, the scattering in the triple-point values of F and 
s are considerably less pronounced. The difference be- 
tween the corresponding values for different interaction 
potentials is comparable to that from different studies of 
a single interaction (see e.g. the data for the Yukawa 
potential). Some (weak) dependence of Stp and Ftp on 
the shape of the interaction potential cannot be com- 
pletely ruled out. At the same time, all the data points 
fall in the relatively narrow range stp — 0.14 ± 0.02 and 
Ftp — 28.0 ± 4.0. This supports the "corresponding 
states principle" conjecture formulated in terms of the 
equality between the values of the parameters s and F, 
characterizing different interacting particle systems. In 
addition, the conditions 

s{p*) ~ 0.14, F(p,,T,)~28, (1) 

are expected to yield reasonable approximations for p* 
and T* at the fcc-bcc-fluid triple point of many other 
interacting particle systems (we will verify this towards 
the end of this Letter). 

In order to illustrate the relation between the corre- 
sponding states principle and the behavior of the in- 
teraction potentials around the mean interparticle sep- 
aration, we have plotted the considered potentials in 
Fig [T] The potentials are normalized by the temper- 
ature, u(x) — U(x)/T, and the reduced distance is 
x = r I A. The parameters of each potential are uniquely 
determined from the conditions s = 0.14 and T = 28 
(which are taken as characteristic for the fcc-bcc-fluid 
triple point). We observe that there is a non- negligible 
difference in the potentials around x = 1, which is ob- 
servable even on the logarithmic scale of Fig [TJ At the 
same time, the first two derivatives are almost coinciding 
in some vicinity around x = 1. 

Next we use numerical results related to the fluid-solid 
and fcc-bcc phase transitions of the considered model 
systems and plot the corresponding phase boundaries in 
Fig. [2] In this figure a portion of the unified phase di- 
agram near the fcc-bcc-fluid triple point is shown in the 
plane of reduced variables s/stp and F/Ftp- The nu- 
merical data are taken from Refs. @, [3 (IPL), [H [H[ 
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FIG. 1: (Color online). IPL, Yukawa, and GCM interaction 
potentials in the vicinity of the mean interparticle separation. 
Potentials are normalized by the temperature, distance is in 
units of the mean interparticle separation, x = r/A. The 
parameters of the potential correspond to s = 0.14 and T — 
28 (see the text for details). The upper (lower) inset shows 
the first (second) derivative of the potentials. 



(Yukawa), and Q (GCM), for details see caption of 
Fig. [5] To evaluate the parameters s and J- along the 
fluid-solid and fcc-bcc phase transitions in IPL and GCM 
systems we used the densities of the fluid and fee phases, 
respectively. For the Yukawa system, the density gap be- 
tween fee and bec phases as well as between fluid and 
solid phases was not resolved in Refs. [22|, [2^1 - The val- 
ues of stp and .Ftp employed are given in Table HI The 
black solid curve corresponds to the universal melting 
equation F(s) ~ 106s 2 / 3 in the regime of sufficiently 
soft interactions (0.1 < s < 1) [1|. 

Figure [2] demonstrates that the phase diagrams of the 
IPL, Yukawa, and GCM systems are topologically equiv- 
alent in the vicinity of the fcc-bcc-fluid triple point. The 
triple point separates the region where the fluid freezes 
directly into the fee crystalline structure (regime of small 
s) from that where the intermediate bee phase is present 
(regime of large s). Another observation from Fig. [5] is 
that the region of stability of the bee phase in the (s, F) 
plane widens when going from the IPL to Yukawa and 
further to GCM interaction. This suggests that the sta- 
bility of the bee phase is promoted by the short-range 
character of the interaction. 

Let us now give several examples of a practical ap- 
plication of the obtained results. Recently, it has been 
suggested that the Lennard- Jones-type 7-6 potential, de- 
fined as U(r) = Ae[{a/r) 7 - (cr/r) 6 ] (where A = 7 7 /6 6 ), 
exhibits freezing into a stable bee structure at sufficiently 
high temperatures [29[. In particular, preliminary results 
from Ref. (29| indicate that the fcc-bcc-fluid triple point 
temperature should be located somewhere in the inter- 
val 5.0 < T* < 100. Since the softness of the LJ-type 
7-6 potential tends to s = 0.143 from below in the high- 
temperature limit, it is very likely that the triple point 
does exist. However, very weak dependence of s on in 
this regime implies extreme sensitivity of the estimated 
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FIG. 2: (Color online). Phase boundaries of the IPL, Yukawa, 
and GCM systems near the fcc-bcc-fluid triple point in the 
plane of reduced parameters s/stp and T/Ttp- Solid sym- 
bols correspond to the fcc-bcc phase transition, open symbols 
mark the fluid-solid (freezing) phase transition. The data 
for the IPL system (circles) are from Ref. Q| (freezing) and 
Refs. 0, [H| (fcc-bcc transition). The data for the Yukawa 
system (triangles) are from Refs. [52] (fluid-solid transition) 
and [24| (fcc-bcc transition). The data for the GCM system 
(stars) are from Ref. @] (both the freezing and fcc-bcc tran- 
sitions). 

triple point parameters to the input value of stp- In 
this case, the criterion ([T]) is not very useful. Taking the 
lower-limit value Stp = 0.12 suggested above, results in 
p* ~ 3.4. From the generic equation relating T* and 
of the n — 6 L J fluids at freezing (Eq. (3) from Ref. |3G |) 
we get T* ~ 32.3. This lower estimate of the triple point 
temperature falls in the interval predicted in [29(. Note 
that the 7-6 L J potential is the only member of the n-6 L J 
family which is expected to form a stable bee structure. 

Another model system displaying interplay between 
solid fee and bec phases is the system of particles in- 
teracting via the exp-6 potential of the form U (r) = 
e {-^=E exp[a(l - r/cr)] - ^ (f) 6 } beyond the hard-core 
at r = ro, such that U(r) exhibits maximum at rrj. We 
concentrate here on the case a = 13, which was consid- 
ered in several studies. Still, the exact location of fcc- 
bcc-fluid triple point is a controversial issue. In Ref. [3l| 
triple point temperature has been estimated as T, ~ 11, 
while in a later study [13 a higher value ~ 20 has been 
reported. Applying the criterion ([!} yields for the triple 
point density and temperature: p* ~ 4.4 and T* ~ 31.3. 
Significant scattering in the triple point temperatures in- 
dicates that further studies are probably warranted to ac- 
curately locate the fcc-bcc-fluid triple point in this model. 

One more potential, which is used to model the re- 
pulsive interactions in systems of soft particles, is the 
so-called Hertzian potential. It is defined as U(r) = 
e(l — r/cr) 5 / 2 for r < a and U(r) = otherwise. This 
potential is bounded (like GCM) and, additionally, has a 
finite interaction range. A detailed study of the phase di- 
agram of Hertzian spheres has been performed by Pamies 
et al. [33d . A complicated picture with multiple re-entrant 
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melting transitions along with first-order phase tran- 
sitions between crystals with different symmetries has 
been reported. Nevertheless, the low density part of this 
phase diagram displays freezing upon compression cither 
into a fee or bec crystal, the fcc-bcc phase transition, 
and the fcc-bcc- fluid triple point. Applying the crite- 
rion ([TJ to the Hertz interaction results in the following 
estimate for the triple point density and temperature: 
p* ~ 1.7 and T* ~ 4.6 x 10 -3 . These values are some- 
what lower than the triple point parameters p„ ~ 1.96 
and T* ~ 7.69 x 10" 3 reported in Ref. [Hj], but still S ive 
an acceptable first guess, which can then be further im- 
proved by detailed numerical simulations. 

Our last example is related to the Yoshida-Kamakura 
(YK) interaction law, U(r) = ecxp[a(l — r/a) — 6(1 — 
r/a) 2 ln(r/er)], where a is a positive parameter control- 
ling the softness of the repulsion. Phase diagram of the 
YK model with a = 2.1 has been reported in Ref (34j 
and with a = 3.3 in Ref. [35| . Again, the reported phase 
diagrams are quite complex, with multiple re-entrant re- 
gions and reach solid polymorphism. The low-density 
part, however, is marked by the presence of the fee and 
bec solid phases, phase transition between them, and the 
fcc-bcc-fluid triple point. The formulated criterion ([T]) 
yields triple-point density and temperature p* ~ 0.24 and 
T* ~ 0.038 for a = 2.1, and p* ~ 0.28 and T* ~ 0.033 
for a = 3.3. These estimates are in fairly good agree- 
ment with the numerical results from Refs. [34| and [35| 
(namely, p» ~ 0.21 and T* ~ 0.040 for a = 2.1 and 
p* ~ 0.25 and T* ~ 0.037 for a = 3.3[36|). 

Thus, the existence and approximate location of the 
fcc-bcc-fluid triple point can be often predicted from the 
knowledge of the first two derivatives of the pair interac- 
tion potential, evaluated at the mean interparticle sepa- 
ration. Although the method is not expected to provide 
very high accuracy, its extreme simplicity can make it a 



useful tool for the preliminary order of magnitude esti- 
mates. In particular, it can be helpful when exploring 
the unique advantages of complex plasmas, associated 
with diversity and variability of interaction mechanisms 
between the macroparticles [fHH, Hi}, in order to visit 
different regions of their phase diagram. Present results 
would be particularly relevant to situations when the 
possibilities to tune and design interactions are present. 
One of the possible methods that has been recently dis- 
cussed [H| is based on the application of external electric 
fields of various polarizations. The physical idea is that 
the applied field induces ion flow, which generates wakes 
downstream from the particles, and thus change inter- 
particle interactions (of course, the frequency should be 
much higher than that characterizing particle dynamics, 
so that the particles do not react to the field). Dif- 
ferent anisotropic interaction classes as well as conven- 
tional isotropic interactions with short-range repulsion 
and long-range attraction can be obtained [381 ] . 

In order to realize these active manipulations of the 
interparticle potential a "PlasmaLab" project has been 
initiated [3Sj . This is a possible followup of the PK-3 
Plus and PK-4 laboratories [4(| for complex plasma 
investigations under microgravity conditions onboard the 
International Space Stations. Two possible plasma cham- 
bers are under active investigations. One chamber has 
a traditional cylindrical geometry with a flexible in- 
ner geometry and electrode system. The other has a 
spherical-like geometry with twelve independently con- 
trolled electrodes. Both chambers are expected to extend 
the plasma parameter space available for investigations 
and support interparticle force manipulations. In this 
context present results could be quite helpful for prelimi- 
nary estimates of the experimental parameter regimes re- 
quired to visit a desired phase state of a complex plasma 
under investigation. 
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